Metallic spin-glasses beyond mean-field theory: An approach to the 
impurity-concentration dependence of the freezing temperature 
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A relation between the freezing temperature (T g ), and the concentration of spins in metallic spin- 
glasses is presented, where T g is obtained as a weighted sum of the absolute values of the spin-spin 
couplings, taking the critical scaling of the correlation function into account. Interestingly, we find 
that no disorder-average is required. We compare our theoretical predictions with experiments and 
find good agreement with the available data from CuMn, AgMn, and AuFe alloys. 

PACS numbers: 



The relationship between the freezing temperature 
(T g ), the distribution of spin-spin couplings [K^), and 
the concentration of impurities (c) in metallic spin-glasses 
has been the subject of vigorous research [lHa. Mean- 
field(MF)-based analytical approaches in ref. [3|,|4j showed 

T g ~ \JHj^[ K Ij] av (where [-] av denotes the average 
over configurations of the disorder). Despite predicting 
the correct low concentration limit T g oc c, this formula 
does not fit the experiments wellQ. A different phe- 

nomenological relation T g ~ J2j=a \J[Kfj]av was derived 
in ref. [H , and an overall good fit to the experiments was 
obtained at moderate c. However, the latter does not 
reproduce the asymptotic linear behavior in the low con- 
centration limit. We derive in this Letter an alternative 
equation for T g that takes the statistical properties of 
the correlations into account, and fits the experimental 
findings well. 

The system: We will be looking at magnetic impuri- 
ties (Fe, Mn), the spins, diluted (0.05 % < c < 7 %) 
in a noble metal (Au, Ag, Cu). At such low concen- 
trations, the spins interact mainly by means of the con- 
duction electrons; i.e., through the RKKY interaction [(| 
K { j cx r~ 3 cos(2fc F r i:( ), where k F is the Fermi wave- 
vector, and is the distance between the spins. The 
quenched positional disorder, together with the strongly 
oscillating sign of K^, results in an even distribution 
of couplings, with zero as dominant mean [HH. Such 
a distribution leads to frustration, and to the formation 
of the spin-glass phase. In this low temperature phase 
(for T < T g ), the spins arc frozen in "random" orienta- 
tions, without a conventional long-range order. In aver- 
age, about half of the spins freezes, with their relative 
orientations trying to satisfy the signs of the pairwise in- 
teractions, whereas the other half is frustrated [l], 0, 0]- 
Two alloys with equal concentration of impurities show 
identical behavior; meaning that the macroscopic state 
is insensitive to the differences among the configurations 
of the microscopic disorder [H, 0, [jj 0] ■ 

The solution: Let us start, considering the spin system 
in the paramagnetic phase, where the thermal average of 



every spin (Sj) is zero. Let S be a test spin, which 
we presume has a non vanishing thermal average (S ) = 
SS . This would change the local fields by Shj = K AQ 8S . 
Shj would affect the neighboring spins, and these, their 
neighbors, and so on. According to the linear-response 
theory, this far-reaching effect is given, at any point, in 
terms of the spin-spin correlations G^ = (S^j), as 5Si = 
j3 J2j Gij 5h^ , with /3 = T~ 1 . Now we relax the constraint 
of Sq having an "externally" fixed thermal average, and 
ask ourselves whether 5S can be maintained by the self- 
induced local fields. If yes, this would imply that the 
paramagnetic phase is unstable, which takes place when 
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More generally, we must consider the effect on all spins, 
starting by variations at an arbitrary number of sites. 
Doing so, we obtain that the paramagnetic phase is un- 
stable when an eigenvalue e Max of Q = f5GK equals one; 
i.e., 

Max(e[Q]) = 1 . (2) 

Eq. {2]) gives the temperature at which the spins ac- 
quire a non vanishing thermal average as a consequence 
of their site-to-site correlations. It holds for spin-glasses, 
as well as for ferromagnets. The specific type of order 
is determined by the actual distribution of K^, which 
also determines the functional form of G^. MF-based 
approaches [H 0, H| , or perturbation expansions Q are al- 
ternatives to avoid the unknown exact form of Gj,-, by 
setting = 5^ at the start, and trying to recover fluc- 
tuations/correlations or introducing reaction fields, at a 
later stage. We proceed here in a different way, and stay 
with Eq.([2]). As the long range of the RKKY interac- 
tion is essential to the glass phase in diluted metallic 
alloys P, [3, [H, a detailed knowledge of the correlations 
between a spin and its nearest neighbors is neither nec- 
essary nor sufficient. We rather look for the general be- 
havior of Gj ■ over the whole interaction range. 

Solving Eq.j2]) seems a formidable problem, but we 
will show now that it leads us back to Eq.{T]). As a 
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consequence of the disorder, neither G,^ nor is a 
self-averaging quantity. The bonds Kj i are uncorrelated 
from each other, and the products G^K-, are weakly cor- 
related. Every matrix element Q u = fiY^j ^ij-^-ji' being 
defined as a sum over all N — 1 sites (j ^ /), samples 
many different K^, covering the whole range of the dis- 
tribution. Then, that Q a is a scalar matrix, and thus self- 
averaging, follows from the theory of large deviations fioj]. 
For i 7^ I, G^ and Kj l correspond to different pairs; the 
product G i jKj l can be either negative or positive with 
equal probability. Then, the sum over j, and therefore 
Q a , converges to zero. For i = I, the sum is greater than 
zero, because G i ■ follows the sign of Kj i as long as the 
frustration allows it, and only about half of the system 
is frustrated [J El, [H, 0]. Qu i s proportional to the lo- 
cal energy density E i = — J2j ^ij^-jf which is position 
independent in a system at thermal equilibrium. In sum- 
mary, all off-diagonal elements are zero, and all diagonal 
elements are equivalent; thus 



(3) 



where i = is an arbitrary origin. The spin-freezing 
occurs, as any other magnetic transition, when the en- 
ergy gained by ordering overcomes the loss of entropy. 
Hence, T g is expected to scale with the interaction en- 
ergy per spin, in accordance with Eq.©. If in Eq.® 
we had done a MF approximation (G^ = 5^ •), we would 
have broken the symmetry of Q a . In such situations, 
authors have resorted to averages over different real- 
izations of the disorder, in order to recover disorder- 
independent equations 0, In contrast to thermal 
averages, the configurational-average does not represent 
a physical process. The actual physical system (with 
quenched positional disorder) does not mutate; it remains 
in one configuration, yet showing sample-independent 
properties. Therefore, there should be an analytical ap- 
proach within which T g can be obtained from a single 
configuration 0, EH; Eq.(j3|) is an example of such ap- 
proach. 

Now we split the sum in Eq.Q into concentric shells 
of thickness dr, with fc F 1 <C dr A; A being the range 
of the interaction, to be defined later. We do first the 
sum in every shell, since this gives us r-dependent aver- 
ages, which will allow us to use the known critical behav- 
ior of the corresponding functions. As half of the bonds 
arc frustrated^, every shell average G^K^ is equal to 
(1/2)16^-1 l-K^I. The positive value, which results from 
the correlations between the sign of (5^5^) and the sign of 
K^, is essential to the phase transition (B), Q- We neglect 
the remaining in-shell correlations between the absolute 
values |G,„| and LfC-,-1, and write 



(4) 



where r c is the typical inter-spin distance, which must be 
proportional to, and not smaller than (3/4-kc) 1 / 3 . Doing 
so, we have not altered the r-dcpcndencc of the integrand, 
and therefore the resulting functional dependence on r c 
and A remains unchanged. Although one would naively 
expect that, in every shell [r;r + dr], larger |GV,|s corre- 
spond to larger |-fTji|s, this is often not the case. Take as 
example any pair of spins for which cos(2k F r i j) = 0. The 
coupling is zero, but G^ is most probably different 
from zero due to indirect correlations. The opposite ex- 
ample is also possible; K ji ^ 0, and G^ = because of 
the frustration. Thus, the effect of the neglected correla- 
tions on the numerical prcfactor is small. However, the 
uncertainty in the value of this prcfactor is unavoidable, 
because the exact distribution of £?_•,■ cannot be calcu- 
lated analytically. 

The critical behavior of\G(r)\: Experiments and Mon- 
tecarlo simulations [l3j agree that the non-linear suscep- 
tibility Xni = (W) Eij (SiSj) 2 diverges at T g , and that 
the spin-glass correlation G SG (r) = (S(r)S(O)) 2 has a 
critical falloff G SG (r) cx (r c /r) 1+n , where i] is the criti- 
cal exponent. As G^ is bounded (i.e., |{S' i S' J )| < 1), the 
double inequality G SG (r) < \G\ < \/G SG (r) holds for 
any r. This implies that \G\ also has an infinite range at 
T = T gl and that 



\G\ « 
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Whether if is a constant (a true critical exponent) or r- 
dependent, is to be found by means of Montecarlo simula- 
tions. A finite cusp in the linear susceptibility at T = T g 
is another characteristic of spin glasses; and it implies 
that the correlation function G has a finite range. The 
long tail of \G\ does not contradict this observation. G 
does not show a critical decay, due to the fluctuating sign 
of Gjj. The double sums J2ij \^ij\ 2 ancl Ylij grow 
faster than the system size (JV), whereas Gij ^ N. 

A simple equation for T g : In order to test Eq.Q 
with available experimental data, we take the scaling 
form ([5]) with rf as a constant, and prefactor 1/2 (uni- 
formly distributed correlations). For \K\, we employ 
a simplified version for the effective RKKY interac- 
tion, due to Shegclski and Geldart 14|: a standard Ir~ 3 
decay with temperature-dependent finite range A T oc 
(TpA/STr/cpT) 1 / 2 > A. A, and T F are, respectively, the 
mean- free-path, and Fermi temperature of the conduc- 
tion electrons. / = 9irJftfi 2 s /T F (2k F ) 3 , where J H is the 
s-d Hund's coupling, and /x efI is the effective moment of 
the impurity. With this, we obtain 



T « 27tc / |G| \K\ r z dr , 



TTCI 

T = 

9 1+7?' 









(- 


' c 




A T 
3 - 


1 



(7) 



3 



Eq.([7]) can be written in terms of directly measurable 
quantities (c, and the resistivity of the sample p), taking 
the form 



TJc = A-B 



4/6 



pTJc 



1+7/ 



(8) 



Eq-© gives the correct behavior T g oc c, in the low 
concentration limit (c — > 0), regardless of the specific 
functional form of p. A and B are given by the follow- 
ing combinations of universal constants and system pa- 
rameters: B/A = [K(3/4 7 r) 1 / 3 (37rne 2 /ftT F ) 1 / 2 ] 1+I '', and 
A = 7iT/(l + 17'). H, and e are the Plank constant and 
the electron charge, respectively; n is the concentration of 
free electrons in the host, k is an universal number, which 
accounts for our uncertainties (i) in the proportionality 
constant in the short length scale r c , and (ii) in the inter- 
In reference 



action range, 
was calculated as 



14|, the effective interaction 
K 2 , instead of \K\. The range of \K\ 

may be slightly smaller than the range of \J K 2 . Thus, we 
should expect r c /A T = K(3/47rc) 1 / 3 (T F A/37rfc F T )- 1 / 2 , 
with 1 < k ~ 10. In the following we will show that 
Eq.® is indeed universal. 

Experimental verification: We start by fitting A, k, 
and rf to the data from the AuFe system, summarized 
by Larsenflij]. This is the only report were c, T_, and 
p{T g ) are tabulated for 17 alloys of a same family (not 
knowing p(T g ), we would have to assume a model for p 
vs c, which would lower the quality of the fit). 
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FIG. 1: The reduced freezing temperature T g /c as a function 
of the damping strength in AuFe alloys. Experimental data 
from Ref.jla] (black dots), and theoretical fit (solid line). 

In FigfTJ we show the experimental data (in black 
dots), together with our best theoretical fit (solid line) 
with A AuFc = 16.8 ± 0.8 K/at.%Fe, B = 2.3 ± 0.5, and 
?/ = —0.41 ± 0.06. Taking the values of the physical pa- 
rameters (T F = 5.5 eV, k F = 1.2 10 10 m -1 , n = 5.9 10 28 
m~ 3 , and p eS ~ 3.25/i B ), we obtain J H w 0.3 eV, and 
k w 10.4. J H falls in the expected range, and it is similar 
to previously reported values derived from the pressure 
dependence of the electrical resistivity 15 1. This indicates 



that our approximation for the prefactor in Eq.® is, al- 
though not very precise, satisfactory. In order to prove 
the universality of Eq.®, and that it really captures the 
physics of the problem, we will verify that we can repro- 
duce the experimental data of other families of alloys, 
employing the values of k and rj' obtained from AuFe. 
To this end, we switch to the AgMn family; data from 
Vier and Schultz (VS)^. 
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FIG. 2: The freezing temperature T g in AgMn alloys vs total 
resistivity, with Sb impurities and 1 at.%Mn (up); vs concen- 
tration of Mn, and no other impurity (down). Experimen- 
tal data from Ref. p^] (black dots), and theoretical prediction 
(solid line). 

The prefactor A AgMn = 6.12 K/at.%Mn is obtained from 
the data point (c = 2.6 at.%Mn, T g = 10.4 K and p = 4 
/if2-cm), where Mn is the only impurity, shown as an 
empty circle in Figj2j Having all three parameters, we 
make a parameter-free prediction of T g as a function of 
p, when non-magnetic Sb-impurities are added (upper 
panel), and as a function of c for clean alloys, taking 
p = p c/c (lower panel). The agreement between pre- 
dictions and experiments is very satisfatory. Relevant 
deviations only appear at high p (A becoming of the or- 
der c -1 / 3 ). We should recall that the expression for the 
effective coupling from Ref. 14| was obtained for a weak 
scattering regime; i.e., low concentration of magnetic and 
non-magnetic impurities. The last two points in the up- 
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per panel correspond to more than 8 at.% of impurities 
(where more than 5% are non-magnetic). 
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FIG. 3: The freezing temperature T g in CuMn allo ys y s con- 
centration of Mn. Experimental data from Ref . [Ia4la | (black 
dots), and theoretical prediction (solid line). 

As a last example, we take the reports on CuMn al- 
loys, from VSpl, (Cowen, Foiles and Shell) [13] , and My- 
dosh et. al.[l8j. The prefactor i CuMn is fitted by forcing 
Eq.© to pass through the point (T g = 9.9 K; p ~ 4 



pQ, cm), of a clean sample CuMn(l at.%) from Ref. [17 1 . 
This point is shown as an empty circle in Fig(3] With the 
obtained value of ^4 CuMn = 16.4 K/at.%Mn, the depen- 
dence of T g on c is predicted (solid line) . Experimental 
data, collected from the literature, are drawn as black 
dots. 

Summary and outlook: I have shown that T g ~ 
(1/2) l^ijll-^jil gi vcs a coherent description of the 

concentration-dependence of the freezing temperature in 
metallic spin-glasses. The sum over sites that naturally 
occurs in this equation, makes unnecessary the aver- 
age over disorder. At the freezing point, |G(r)| has a 
long tail, cx r~^'\ with -1 < 2=i < rf < rj. The 
fit to the experimental data gives if ps —0.41. This 
seems to rule out the hypothesis of a chiral glass tran- 
sition (77 ps 0.6) [Hj], leaving open the choice between 
the Ising glass (77 = —0.41) and the Heisenberg glass 
(77 = —0.3) alternatives. Studying |G(r)| by means of 
Montecarlo simulations and performing more measure- 
ments of (c; T g ] p(T g )) would be of great help in finding 
the exact value (s) of 77'. It may be also interesting to 
look for possible cross-overs (changes in 77') as a func- 
tion of the concentration and type of impurities, specially 
when mixing elements with different spin-orbit coupling 
strengths. 

The author thanks professor Mark R. A. Shegelski for 
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